%setsubset

%all = {[1:53]};
all = [1:32 34:53];
real = [2:6 9:11 14:15 17 19 21 23:32 34:36 38:39 41 45:46 48:49 51:52];
%ALL_COAL_1 = {[1:9 11:53]};
%    AGR = {[1 12 14:21 33 38:41 45:53]};
    agr = [1 12 14:21 38:41 45:53];
    nonagr = [4:11 13 22:32 34:36 42:44];
%REAL = {[2:6 9:11 14:15 17 19 21 23:36 38:39 41 45:46 48:49 51:52]};
 %       AGR_REAL =  {[14:15 17 19 21 33 38:39 41 45:46 48:49 51:52]};
        agr_real = [14:15 17 19 21 38:39 41 45:46 48:49 51:52];
        nonagr_real = [4:6 9 11 23:32 34:36];
            nonagr_real_phys = [22:32 34:35];
            nonagr_real_nonphys = [4:6 9 11];
    nom = [1 7 8 12 13 16 18 20 22 37 40 42:44 47 50 53];
        agr_nom = [1 12 16 18 20 37 40 47 50 53];
        nonagr_nom = [7:8 13 22 37 42:44];

